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ABSTRACT 

m ; 

Excess reddened optical emission has been measured in the power-law decay 
afterglow light curves of several GRBs some tens of days after the GRB events. 
O ' This emission is claimed to be a signature of a supernova (SN) taking place at 

the same time as the GRB event, in support of the collapsar scenario. Obser- 
t— i ■ vations of features in prompt and afterglow X-ray spectra of GRBs require high 

' density, high metallicity material near GRB sources, supporting the supranova 

(SA) scenario. These conflicting lines of evidence are resolved in likely regimes of 

■ parameter space if the supernova remnant (SNR) shell is illuminated and heated 
by a pulsar wind (PW), and cools following the GRB event. The heating of 
the SNR shell by the PW is parameterized, and simplified expressions for the 
synchrotron, synchrotron self-Compton, and thermal radiation fields within the 
SNR shell and pulsar wind bubble are derived for monoenergetic injection of wind 

J^j ■ electrons. A cooling SNR shell could produce the excess emission detected from 

GRB 970228, GRB 980326, and GRB 011121. X-ray synchrotron radiation from 
the PW electrons provides an important source of ionizing radiation. 

Subject headings: gamma-rays: bursts — gamma-rays: theory — radiation pro- 
cesses: nonthermal — pulsars 



1. Introduction 

Two leading scenarios to explain the origin of gamma-ray bursts are the collapsar and 
supranova (SA) models (see Meszaros (2002); Dermer (2002) for review). The collapsar 
model (e.g., Wang & Woosley 2002; Woosley 1993) assumes that GRBs originate from the 
collapse of a massive star to a black hole. During the collapse process, a nuclear-density, 
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several Solar-mass accretion disk forms and accretes at the rate of ~ 0.1-1 M s _1 to drive 
a baryon-dilute, relativistic outflow through the surrounding stellar envelope. The duration 
of the accretion episode corresponds to the prompt gamma-ray luminous phase, which is 
commonly thought to involve internal shocks. A wide variety of collapsar models can be 
envisaged (Fryer, Woosley, & Hartmann 1999), but their central feature is the one-step 
collapse of the core of a massive star to a black hole. 

Optical observations of GRB afterglows reveal significant evidence for enhanced emission 
in excess of an extrapolation of the optical power-law decay spectra in a few cases, namely 
GRB 980326 (Bloom et al. 1999), GRB 970228 (Reichart 1999; Galama et al. 2000), and 
GRB 011121 (Bloom et al. 2002; Price et al. 2002). The spectrum of the excess is reddened 
at frequencies above the spectral peak, as might be expected for a thermal emitter. The 
peak of the excess emission occurs tens of days after the GRB event (in the local frame), 
and the SN event precedes GRB 011121 by ~ 3-5 (±5) days (Bloom et al. 2002), consistent 
with a one-step collapsar model. 

The supranova model (Vietri & Stella 1998) of GRBs involves a two-step collapse process 
of an evolved massive star to a black hole through the intermediate formation of a neutron 
star with mass exceeding several Solar masses. The neutron star is initially stabilized against 
collapse by rotation, but the loss of angular momentum support through magnetic dipole and 
gravitational radiation leads to collapse of the neutron star to a black hole after some weeks 
to years. The accretion-induced collapse of a neutron star in a binary system could also 
form a GRB (Vietri & Stella 1999). A two-step collapse process means that the neutron star 
is surrounded by a SNR shell of enriched material at distances of ~ 10 15 -10 17 cm from the 
central source. The earlier SN could yield ~ O.1-1M of Fe in the surrounding vicinity. The 
discovery of variable Fe absorption during the prompt emission phase of GRB 990705 (Amati 
et al. 2000) and X-ray emission features in the afterglow spectra of GRB 991216 (Piro et al. 
2000) has been interpreted in terms of the SA model (Lazzati et al. 2001; Bottcher et al. 2002; 
Ballantyne et al. 2002). Low significance X-ray absorption features observed in GRB 970508 
(Piro et al. 1999), GRB 970828 (Yoshida et al. 2001), and GRB 000214 (Antonelli et al. 
2000) can also be modeled within the SA scenario, preferably with photoionization reflection 
spectra (Lazzati et al. 1999). The multiple features detected in GRB 011211 (Reeves et al. 
2002) are also suggested to arise from an excited SNR shell at ~ 10 15 cm from the central 
source. 

A PW and pulsar wind bubble (PWB) consisting of a quasi-uniform, low density, highly 
magnetized pair-enriched medium within the SNR shell is formed by a highly magnetized 
neutron star during the period of activity preceding its collapse to a black hole (Konigl 
& Granot 2002). The interaction of the PW with the shell material will fragment and 
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accelerate the SNR shell, and the PW emission will be a source of ambient radiation that 
can be Comptonized to gamma-ray energies (Inoue et al. 2002; Guetta & Granot 2002). 

Here we propose that the nonthermal leptonic PW radiation provides a source of hard 
ionizing radiation to heat and photoionize material in the SNR shell, which is described in 
terms of a partial covering model as a result of the strong PW/SNR interaction. This radia- 
tion heats the shell, and the escaping thermalized radiation is seen as the delayed reddened 
optical excesses in GRB optical afterglow light curves. The source of ionizing nonthermal ra- 
diation, which is extinguished when the GRB occurs, would provide a nonthermal precursor 
to a GRB, though at low flux. Various other features of this system are outlined. 

2. Pulsar Wind Physics 

A very rapidly rotating, strongly magnetized neutron star is formed in the first step of 
the SA scenario. For a neutron star with surface polar field of 10 12 Bi 2 G, a rotation rate of 
10 4 f24 rad s _1 , and a radius of 15ri 5 km, the radiated power L = 7 x 10 44 B 2 2 r\ b Q,\ ergs s _1 . 
The quantity E rot = \j(GM 2 /c 2 )tt = 2 x 10 53 jo.sMf^ ergs = 10 53 £^' is the initial 
rotational energy of the neutron star, expressed in terms of the neutron-star mass M = 
3M 3 M and the dimensionless angular momentum j .5 — GM 2 /2c, which remains roughly 
constant during collapse (Vietri & Stella 1998). The spin-down time t sd = (E rot /L ) = 
2.8 x lO^o.sMK^VM)- 1 s. 

The torque equation, (l = —K£l n , where the braking index n = tlil/(l 2 , can be solved to 
give Q(t) = Q (l + t/td)~ 1 ^ n ~ 1 \ provided that K remains constant (e.g., Atoyan 1999). (For 
the Crab, n = 2.5.) We assume that the spin-down power is radiated as a PW in the form 
of particle and Poynting flux, so that the wind luminosity L w = L sd = (£ e + £ p + ^b)L S( i = 
L /(l + t/t s d) k , where k = (n+ l)/(n— 1), and we have divided the wind power into leptonic 
("e"), hadronic ("p"), and electromagnetic ("B") fractions £j, % = e,p,B. 

In the SA model, rotational stability is lost due to spindown, and the neutron star 
collapses to a black hole to produce a GRB at time t co u- Here we assume t co u < t s d, so that 
for k > 1 the wind energy E w (t) = i] (k - l) _1 L t S(i [1 - (1 + (t/t sd ) l - k ] -> r]L t when t < t sd , 
where the parameter 7] roughly accounts for the escape of wind energy from the SNR shell 
as well as losses of wind energy due to radiation, and is assumed here to be constant with 
time. The equation for the SNR shell dynamics in a spherical approximation for the SNR 
shell is MsnrR = 47iR 2 p w (t), where Msnr = it^snrMq is the SNR shell mass, the wind 
pressure p w (t) ~ E w (t)/3V pw b, the PWB volume V pw b = 47ri? 3 (t)/3, and the PWB radius 
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R(t) = J Q dt' v(t') (Inoue et al. 2002). This equation can be solved to give 



R(t) = „ of (l + JL) , where U = ~ 1.9 X 10' ^» s (1) 

(Guetta & Granot 2002). It follows that t acc /t sd = E^ R /r]E rot ~ 0.lE^ NR /r]E^, where 
^Msjvflug = 9 x 10 51 m S jvfl/3ii = i0 52 E§f R ergs is the SNR shell kinetic energy, and w = 
0.1/3_ic is the SNR shell coasting speed. If 77 = 1, then t acc <C t s a for standard values, 
and shell acceleration must be considered. A highly porous shell has rj <C 1. This will 
occur for strong clumping of the ejecta, which takes place on the Rayleigh- Taylor timescale 

t RT ~ y/ x/R, where the characteristic clumping size scale is x = fR (Inoue et al. 2002). 

Hence (tux /tacc) = f 1 ^ 2 (t/t acc ) 3 ^ 4 in the regime t < t acc . This shows that small-scale (/ <ti 
R) clumping will occur when / \/t/t acc . Only a detailed hydrodynamic simulation can 
characterize the porosity of the shell due to effects of the pulsar wind, which 77 parameterizes. 
Here we assume, in contrast to the picture of Konigl & Granot (2002), that a large fraction 
of the wind energy escapes the SNR shell and rj ~ 0.1, so that shell acceleration can be 
neglected. This effect also causes the shell to become effectively Thomson thin much earlier 
than estimated on the basis of a uniform shell approximation, except for high density clouds 
with small covering factor. A general treatment is needed to treat shell dynamics, including 
travel-time delays of the pulsar wind and deceleration of the SNR shell by the surrounding 
medium. When t co u < t acc , t s d, as assumed here, R = v t. 



3. Radiation Field in SNR Cavity 

The strong PW provides a source of nonthermal leptons, hadrons, and electromagnetic 
field. Dominant radiation components considered here are leptonic synchrotron and SSC 
radiation, and thermal emission from the inner surface of the SNR shell which is heated by 
the nonthermal radiation. 



3.1. Wind Synchrotron Radiation 

The volume-averaged mean magnetic field B w in the SNR cavity powered by the PW is 
obtained by relating the wind magnetic field energy density ub = B^(t)/8ir = eB W u w (t) = 
?>eB W ijLQt/i.T{R? through the magnetic- field parameter e^, also assumed to be constant in 
time. Thus 

5«(f) = f, where B = = 4 x 10 8 v /^ (^£1) G-s. (2) 
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Let 7 W = 10 5 75 represent the typical Lorentz factors of leptons in the wind zone of the 
pulsar. This value should remain roughly constant when t < t sd . A quasi-monoenergetic 
proton/ion wind may also be formed, though we only consider leptonic processes here, and 
furthermore do not treat nonthermal particle acceleration at the PW/SNR shell boundary 
shock. When synchrotron losses dominate, the mean Lorentz factor of a distribution of 
leptons with random pitch angles evolves in response to a randomly oriented magnetic field 
of mean strength B according to the expression —^ syn = [[a^B 2 / 'Qimi e c)^ 2 , giving 

1 W = [\-T B {\-\)Y 1 , (3) 

7 Li L 

where ti is the injection time, 7, = 7^ is the injection Lorentz factor, and 

T Blw = 7 ^2x 10 13 e BwV (^^) 75 s . (4) 

Writing the nonthermal electron injection function as 

dNej^U) = ri( e L _ ^ 
dUd^fi m e c 2 ^ w 
the electron Lorentz factor distribution dN e (j] t)/d r ) can be solved to give 

^v&O = ^ 1 1 _ 

07 t 7 

In this expression, Af e = r)£ e L / '(m e c 2 ^ w ), and we introduce the dimensionless quantities 
t = t/Ts^w and 7 = 7/7^. Adiabatic losses, which are significant on time scales ~ t/3, 
are negligible in comparison with synchrotron losses of wind electrons when t <C 1, and we 
restrict ourselves to this regime. 



The vL v synchrotron radiation flux vL s ^ n ~ \u b cgt?^ N e ("y s ), where 7 S = \fefe~B, 
e = hv/m e c 2 , ts = B/B cr , and the critical magnetic field B cr = 4.41 x 10 13 G. Thus 



vLT{e) = car ^ (KT Bl 2 J {jf 2 ^(l + J^f- 1)~ 2 , for e < e max , (7) 
l07r t t V et 



where u = B/ lw B cr T B = 4.5 x 10~ 19 p^ /(^Tje^B 12 r 3 15 Qh 5 ), and 

B_ 4.4 x IP' 9 /g^Ts 



e 6 7^, 4.4X10-° 



The synchrotron cooling timescale for wind electrons is t syn = ^wTsi 2 . In the regime 
t < Po where the wind electrons strongly cool, we can approximate the synchrotron spectrum 
by 

vLT\e) = h e L (— ) 1/2 H(e; e , e max ) , (9) 
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where H(x; a, b) = 1 for a < x < b, and H = otherwise. The value of i at t s d is 

tsd ~i4x io- 5 ^ M 3 2 l-i (10 ) 



Strong cooling holds when t < t s d/T B ^ w1 noting that time in physical units can be found 
from equation (4). 



3.2. Wind Synchrotron Self-Compton Radiation 

The relative importance of synchrotron self-Compton (SSC) to synchrotron cooling is 
given by the ratio p = u syn /uB of the synchrotron radiation energy density to Mg = B 2 /8nt 2 
(see eq.[2]). The PW synchrotron radiation energy density u syn ~ K(vL™ ax )/4:7rR 2 c, where 
v jmax ^ £ e L and the parameter k > 1 is a reflection factor such that k = 1 corresponds to 
perfect absorption or direct escape (small covering factor) of the PW synchrotron radiation, 
and k 3> 1 corresponds to a highly reflecting SNR shell with large covering factor. One 
obtains p = 0.033/t£ e/ 3_i/e£«)f7. For large porosity with k = 1 and r\ pa 0.1, we see that 
the SSC contribution is small compared to the synchrotron flux when i e f3-i/e Bw ^ 1- We 
restrict ourselves to this regime where equation (3) is valid. 

In the 5-function approximation for the energy gained by a photon upon being scattered 
by relativistic electrons (Dermer, Sturner, and Schlickeiser 1997), with scattering restricted 
to the Thomson regime, the vL v SSC spectrum from the cooling wind electrons is given by 

v - Snv 2 T B i 2 Vw [ a (ao + x) +a ° [H a + x )l ^ W 

where ao = i^ 1 — 1, xq = 7«,\/eo/e, and x\ = 7 W A/max(e, e^ 1 , e max )/e, where we have used 
approximation (9) for the synchrotron spectrum. 



3.3. Thermal Emission from the Interior of the Shell 

The SNR shell is pictured as having been shredded by the PW and therefore highly 
porous, though permeated with small scale (~ 10 12 -10 14 cm) density irregularities consisting 
largely of metals (Dermer & Mitman 1999; Bottcher et al. 2002). Note that O.lm 56 _iM 
of 56 Ni provides a radioactive decay power of ~ 2 x 10 42 to 56i -i ergs s -1 before decaying on 
an ^ 120 day timescale. The partial covering factor p c (< 1) corresponds to the covering 
fraction by optically-thick, dense shell inhomogeneities. Recognizing the severe limitations 
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in the following expression, we determine the effective temperature T e ff in the interior of 
the shell through the relation 

47r J RV SB T e 4 // p c A Pc ^ + 2 x 10 42 m 56 ,-i exp(-t/120 d) ergs s" 1 , (12) 

where *4.(< 1) is an absorption coefficient, dependent in general on the evolving radiation 
spectrum, though A is here assumed constant. The wind heating is more important than 
heating from the decay of radioactive Ni when (A/ 0.1) (p c / 0.1) (£ e /l/3) > m 56 ,_i/ 1 Bl 2 r\^Vl\. 
In this regime, 

^ W SB j - ^ 4/3 ) V^^M 3 K ' (13) 

letting k s; 1 because p c ~ 0.1. For standard parameters, the shell is heated to ~ 10 4 K 
temperatures when t/t s d <C 1. When B i2 > 4, this occurs at t/t sd < 1. The PW heat 
source is extinguished once the GRB occurs, and the thermal emission from the SNR shell 
will decay on the dynamical timescale t dyn = R(t)/c = 324(3_i (j ^ 5 M^ / Bf 2 rf 5 Ql) (t/t sd ) d, 
noting that the cooling timescale of a hot shell is generally shorter than t dyn . 



4. Results and Discussion 

Fig. 1 shows the total vL v flux composed of synchrotron, SSC, and thermal components, 
calculated according to the above simple analysis. Standard parameters are j _ 5 = B 12 = 
r 15 = fi 4 = (3_i_ = k = M 3 = m 52 „i = 75 = 1, £ e = e Bw = 1/3, and p c = r] = A = 0.1. 
The heavy solid and dashed curves correspond to standard parameters with t/t sd = 0.1 and 
0.01 as labeled, with t dyn = 32.4 and 3.2 days and T e jf = 1700 and 5420 K, respectively 
(t sd = 3240 d). This system is in accord with the standard SA model and might explain 
the delayed excess emission in GRB 011121, contrary to the conclusions of Bloom et al. 
(2002). The set of three dotted curves employ standard parameters, except that 5 12 = 4, 
and t/tsd = 1, 0.1, and 0.01 as labeled (t sd = 202 d), with t dyn = 20.2, 2.0, and 0.2 days (= 17 
ks), and T eff = 3760, 11900, and 37630 K, respectively. 

The intense PW X-ray synchrotron emission provides a source of ionizing photons 
throughout the SNR shell that would decay on the timescale t dyn (the line flux would decay 
on si 2t dyn ). The B u = 4, t/t sd = 0.01 case might correspond to a GRB 011211-type system. 
The PW radiation field provides an additional source of ionizing photons in the SA model 
not considered by, e.g., Kumar & Narayan (2002), and would be energetically important for 
this case if 75 = 0.1. A varying ionization flux impinging on an expanding SNR shell will 
produce a characteristic kinematic signature of shell illumination. Our results differ from 
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the conclusions of Guetta & Granot (2002), because we assume that the SNR shell highly 
porous, as in the scenario of Inoue et al. (2002). 

The thin solid curve is the result for standard parameters except that B u = 2, 7^ = 10 6 , 
and t = 0.1t sd = 81 d. For this case, t dyn = 8.1 d and T e ff = 4360 K. It is apparent that 
a sizeable parameter regime could make quasi-thermal emission between ~ 10 14 -10 15 Hz at 
the level of ~ 10 42 -10 43 ergs s _1 some tens of days after the GRB event. 

Thus, even if the GRB takes place months after the SN, a delayed reddened excess from 
the cooling shell may be seen, for plausible parameters, some tens of days after the GRB. 
The PW heating is in addition to the the heating due to the Ni— > Co — > Fe chain, and may 
relax the need to have ^56,-1 ^> 1, as required to explain the unusually bright light curve of 
SN 1998bw (Iwamoto et al. 1998). The thermal emission from the SNR shell heated by the 
PW is proposed to be the origin of the reddened optical excess emission detected in GRBs 
970228, 980326, and 011121. 

This model predicts decaying nonthermal ionizing hard X-ray emission, which could be 
detected with INTEGRAL or Swift from nearby z < 0.1-0.3 GRBs when the PW synchrotron 
flux is at the level > 10 ergs s^ 1 . A nonthermal optical emission signature would accompany 
the reddened excess, with a relative intensity that depends sensitively on 7„. A hard X- 
ray sky survey instrument such as EXIST could monitor for the PW nonthermal emission 
preceding aligned and off-axis GRBs. The X-ray synchrotron radiation provides a source of 
external radiation to enhance photomeson production by energetic hadrons over the standard 
GRB model (Dermer 2002; Atoyan & Dermer 2001). High-energy neutrino and GLAST 
gamma-ray observations can distinguish the collapsar model and different versions of the SA 
model (A. M. Atoyan & C. D. Dermer, 2002, in preparation). 

I would like to thank M. Bottcher, A. M. Atoyan, and S. D. Wick for discussions and 
comments. This work is supported by the Office of Naval Research. 
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Fig. 1. — Total model spectral energy distribution composed of synchrotron, SSC, and shell 
thermal emission from a PW and wind-heated SN shell. Standard parameters, given in the 
text, are used, except where noted in the legends. 



